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ABSTRACT 

> 

^ \ We model SNRs at a variety of heights above the disk with a detailed nu- 

merical simulation that includes non-equilibrium ionization and recombination 
and follows the remnants' evolution until their hot bubbles have cooled. We an- 
^ \ alytically calculate the bubbles' buoyant acceleration and frictional drag. From 

the simulation results, combined with the rates for isolated supernova explosions 
above a height of 130 pc, we estimate the time and space average O"*"^, N"*"^, and 
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S I C"*''^ column densities and emission intensities, 1/4 keV soft X-ray surface bright- 



O ■ 

^ \ ness, area coverage, and volume occupation due to the population of isolated 

SNRs above the Galaxy's HI layer. Irrespective of assumed supernova explosion 
energy, ambient nonthermal pressure, or frictional drag coefficient used in the 



X 



\ calculations, the predicted O"*"^ column density as a function of height matches 

the observed distribution between 130 pc and 2000 pc. The O VI resonance line 
emission (1032, 1038 AA) contributes significantly to the average observed inten- 
sity. Assuming our modest supernova explosion rate, the population of isolated 
extraplanar SNRs can explain 80% of the observed 1/4 keV surface brightness at- 
tributed to the extraplanar gas beyond the H I layer in the southern hemisphere. 
Within the range of uncertainty in the SN rate, such SNRs can explain all of this 
observed emission (400 counts s~^ arcmin"^). Thus, extraplanar SNRs could be 
the most important sources of hot gas between the Local Bubble and z ~ 2000 pc 
in the relatively quiescent southern hemisphere. These results stand whether the 
remnants are assumed to be buoyant or not. The population of old extraplanar 
SNRs should cover most, but not all of the high latitude sky, thus explaining 
the mottled appearance of the soft X-ray maps (outside of superbubbles) . Bright 
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young extraplanar SNRs should cover less than 1% of the high latitude sky. Per- 
haps the I — 247°, b — —64° crescent in the 1/4 keV X-ray maps could be such a 
remnant. 

Subject headings: hydrodynamics - Galaxy:general - ISM:general - ISM:supernova 
remnants ultraviolet: ISM - ultraviolet: ISM 



1. Introduction 

About half of the progenitor stars at the Sun's Galactocentric radius reside above the 
Galactic H I disk (Ferriere 1998). After they explode, their remnant bubbles will evolve 
in a relatively rarefied, dust-poor environment. They will expand to greater sizes, have 
lower interior pressures, lower luminosities, lose less energy via emission by dust grains, and 
live longer than disk remnants (Cioffi 1991). At any given time, most of the supernova 
remnants (SNRs) above the H I layer should be dynamically old, having transitioned from 
the adiabatic to the radiative phase in the first few percent of their lifetimes. Old remnants 
should contribute to the 0+^ and soft X-ray backgrounds even though they are dim (Shelton 
(1998, 1999), hereafter Papers I and II). The young remnants should be much brighter 
and so easier to recognize. Known examples of comparatively young extraplanar supernova 
remnants include SN1006 at 2; ~ 450 pc (Laming et al. 1996) and the Lupus Loop at 
z ~ 330 pc (Leahy et al. 1991). In addition, bright arcs in the ROSAT 1/4 keV maps at 
/ = 247, b = —64 and / = 215, b = —68 may also be extraplanar SNRs. 

The plenitude of isolated SNRs above the Galactic H I layer raises the possibility that 
they play important roles in heating and ionizing the extraplanar interstellar medium (ISM) 
and whets our appetite for more information. Considering that SNRs evolving in low density 
environments can be longlived, might they cover a large solid angle on the high latitude sky 
or fill a large fraction of the volume of space? Might they create many of the observed high 
latitude high-stage ions (O"^^, N+^, and C+^) and soft X-rays? Considering that young SNRs 
are both bright in soft X-rays and rich in O"^^ ions, do soft X-rays and 0+^ ions track each 
other in general, or only in recently shocked gas? How buoyant is the hot gas within SNRs? 

Much progress has been made toward understanding these SNRs. According to Cioffi 
(1991) 's and Paper I's rough estimates, SNRs fill a small fraction of the extraplanar volume, 
but cover a large fraction of the high latitude sky. ^ Paper I's preliminary estimate for the 



^This paradox is possible when a large number of SNRs exist at any given time and each is small on 
the scale of the available space (fractional area coverage goes as N{R/L)'^, but fractional volume occupation 
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number of O"^^ ions produced by the population of isolated, extraplanar SNRs accounted 
for the extraplanar O^^ observed by Copernicus (Jenkins 1978a) and the supernova explo- 
sion scalchcight (300 pc) matched the observationally determined 0+^ scaleheight (Jenkins 
1978b). Subsequently, Ff/^i? observed a much larger sample of high z stars and found more 
O^^ above 300 pc. These new observations suggest that modeling the vertical distribution 
of SN explosions may be important. FUSE also found O"*"^ moving away from the Galactic 
plane at more than 100 km sec~^ (Wakker et al. 2003). This suggests that buoyancy may 
be important. Another possible sign of buoyant motion in the interstellar medium is an 
H I "mushroom-shaped cloud", whose cap may be the squashed shell of a buoyant SNR 
and whose stem may consist of traihng material (English, et al. 2000). Multidimensional 
simulations have been able to reproduce mushroom shapes and buoyant acceleration (Jones 
1973; Avillez & Mac Low 2001). Could buoyancy explain the large 0+^ column densities 
and scale height? 

If we could envision a map of the sky based on the O^^ column density in extraplanar 
SNRs, it would not look like a similarly constructed map of O"*"^ resonance line emission, 
or, for that matter, like a similarly constructed map of soft X-ray emission. The O"*"^ ions 
would be more smoothly distributed and cover greater area than the detectable emission. 
The reason is that SNRs harbor rich stores of high-stage ions throughout their lives, but are 
bright emitters only during their early stages when the dense, recently shocked gas behind 
their shockfronts is hot (Slavin & Cox 1992; Cioffi 1991; Shelton 1998, 1999). Probably 
other types of hot gas structures behave similarly. This difference may explain the difficulty 
in matching the FUSE 0+^ column density map (Savage et al. 2003) with the ROSAT 1/4 
keV map (Snowden et al. 1997). Furthermore, extraplanar SNRs are better emitters of 1/4 
keV photons than 3/4 keV photons. Thus, the more smoothly distributed 3/4 keV emission 
(McCammon ct al. 1983; Snowden ct al. 1997) cannot be attributed to isolated SNRs, even 
though some or all of the lumpy 1/4 keV emission can (Paper 11, Slavin et al. (2000)). 

The Paper I and II analyses were preliminary in the sense that they relied on simulations 
of SNRs born at a single height above the plane and were meant to be followed by a suite 
of SNR simulations. The Slavin et al. (2000) project used the entire SNR population rather 
than the population of isolated SNRs above the Milky Way's thick disk. Furthermore, 
their simulations did not include ambient magnetic pressure which can compress old SNRs 
and therefore increase their UV and X-ray luminosities. Consequently, these two studies 
predicted significantly different time and space average 1/4 keV surface brightnesses. We 
resolve these difficulties in the present paper. 



goes as N{R/LY, where N is the number of objects and R/L is the object's radius divided by the length 
scale of the available space, thus R/L < 1.) 
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Here, we aim to estimate the observable effects of the population of extraplanar SNRs 
by simulating a suite of remnants at various heights above the Galactic plane and various 
magnetic field strengths and explosion energies, allowing for buoyancy, and by multiplying by 
the expected rate of isolated supernova explosions above the Galactic disk. The simulation 
package used to calculate the time evolution of isolated supernova remnants is described in 
detail in Paper I and described briefly in Subsection 2.1. The analytic calculations used to 
predict the buoyant acceleration are spelled out in detail in Subsection 2.2. The ambient 
conditions and progenitor statistics are presented in Subsections 2.3 and 2.4. Due to the 
interstellar medium's variation with height above the plane and due to the uncertainty in 
the explosion energy and ambient non-thermal pressure, we ran a suite of simulations. Their 
parameters arc presented in Subsection 2.5. The resulting area coverage, volume occupation, 
Q+5^ N"*"^, C+'^, and soft X-ray predictions for each of the simulated SNRs are tabulated in 
Subsection 3.1. We combined the individual SNR predictions with the buoyancy calculations 
and convolved with the rate for isolated supernova explosions above the Galactic disk's 
H I layer. The resulting estimates of the time and space averaged area coverage, volume 
occupation, high stage ion content and emission, and 1/4 keV X-ray intensity due to the 
population of extraplanar SNRs are presented and discussed in Subsections 3.2 - 3.5. In 
Section 4, we compare the estimates with O"*"^ and 1/4 keV X-ray observations. O"*"^ and 
soft X-ray emission are far less affected by photoionization than are C"*"^ and N"*"^, so are 
better tests on the model. The column densities from the calculations agree well with the O"*"^ 
column densities observed within the first 2000 pc of the Galactic plane. The simulations also 
explain much of the observed extraplanar 1/4 keV X-ray emission outside of superbubbles 
and explain its mottled appearance. We summarize the results in Section 5. 



2. Simulation Method 

2.1. Computer Simulation Package 

The supernova simulations are performed with a Lagrangian mesh hydrocode whose 
algorithms model shock dynamics, nonequilibrium ionization and recombination, and non- 
thermal pressure. Previous investigations found that some form of mixing is necessary in 
order for the remnants' centers to retain the modest gas densities surmised from X-ray obser- 
vations of post Sedov phase SNRs such as W44. The hydrocode models thermal conduction, 
which diffuses entropy within the remnant's interior and so approximates the required mix- 
ing (Cox et al. (1999); Shehon et al. (1999); Shelton, Kuntz & Petre (2004)). Cui & Cox 
(1992) found that the electrons and ions behind the SNR shock zone come into equilibrium 
early in the remnant's evolution, long before radiative cooling becomes important. As a re- 
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suit, the SNR's late term evolution (which plays the largest role in calculating the remnant's 
O"*"^ content, area coverage, and volume occupation) is not affected by our approximation 
of instantaneous electrons and ion equilibration. Cosmic rays affect SNRs by contributing 
to the pressure behind the shock, while SNRs affect cosmic rays by accelerating them. Cos- 
mic rays affect SNRs by contributing to the pressure behind the shoclc, while SNRs affect 
cosmic rays by accelerating them. In a careful analysis of supernova remnant evolution in 
the presence of cosmic rays and magnetic fields, Ferriere & Zweibel (1991) explored several 
detailed interactions, but, nonetheless, found that their simulated SNRs with cosmic rays 
were qualitatively similar to their simulated SNRs without cosmic rays. In our simulation 
package, we include an isotropic cosmic ray pressure (see Section 2.5), but neglect detailed 
interactions. Furthermore, we assume that the interstellar medium's gas phase metal abun- 
dances are approximately solar and so adopt the abundances of Grcvcssc & Anders (1989). 
The code is similar to Slavin & Cox (1992, 1993)'s code. Papers I and II used this code, 
describe it in detail, and present sample simulations. Because the hydrocode's Lagrangian 
mesh architecture does not model a stratified ambient media, we combine the SNR simu- 
lations with the following analytic buoyancy calculations in order to determine if and how 
rapidly the remnants rise in the Galaxy's stratified interstellar medium. 



2.2. Buoyancy 

As a whole, a SNR is not buoyant, because its average density approximately equals that 
of the surrounding gas. Only the very hot, rarefied interior should be buoyant. Furthermore, 
buoyancy is probably not relevant for very young remnants because the buoyant acceleration 
acting on the bubble interiors is outpaced by the remnants' expansion and because their 
magnetic fields tie the most buoyant material in their centers to the non-buoyant dense, 
ionized, postshock gas. Later in a remnant's evolution, a cool shell develops between its 
hot bubble interior and its shockfront (see Slavin & Cox (1992). or Paper I for detailed 
simulations). At this time, buoyant motion may exceed the hot bubble's expansion and the 
postshock shell cools, becomes somewhat neutral, and dissipates. It is in this stage that 
we consider the buoyant force. Because we do not know how well the shell constrains the 
hot interior nor the strength of frictional drag, we consider a range of cases. In the first, 
we assume that the SNRs do not rise buoyantly. In the second, we assume that after the 
cool shell develops, the hot remnant interiors dissociate from the shell and rise buoyantly, 
but encounter some resistance. In the third, we assume that after the cool shell develops, 
the hot remnant interiors dissociate from the shell and rise buoyantly, without resistance. 
Below, we describe our analytical calculations for the buoyant motion of a hot bubble interior 
encountering frictional drag. 
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The elementary equation for the acceleration due to buoyancy is: 



Ph[z) 

where g is the acceleration of gravity, pa is the mass density of the ambient medium, pf, 
is the mass density within the bubble, and z is the height above the Galactic plane. The 
acceleration of gravity at the Sun's galactocentric radius is given by Benjamin & Danly 
(1997) and Wolfire et al. (1995) as: 

g{z) = -9.5 X 10"^ tanh(^/400pc) cm s'^ (2) 



As the bubble gains speed, it should experience increasing frictional drag. Following 
Benjamin & Danly (1997), we calculate the frictional acceleration in the direction opposite 
to motion from ^ 

a/ = -CAbv'^pa/mb, (3) 

where is the cross sectional area of the hot bubble, Vb is the hot bubble's velocity relative 
to the ambient medium pa is the ambient mass density, rub is the mass of the hot bubble 
(determined from following simulations) and C is the drag coefficient. If the bubble slides 
frictionlessly through the ISM, then C — 0. However, if the bubble sweeps up the overlying 
ISM, then C — 2. If, in addition, a low pressure region develops behind the bubble, then 
C exceeds 2. The value of C for hot bubbles in the ISM is poorly known, but simulations 
suggest (Jones et al. 1996) that it may be around 1. Here we consider 3 cases, C — 0, C — 1, 
and C — large enough to prevent bubble motion. 

Generally, the center of the hot bubble is the most rarefied and thus most buoyant, while 
the periphery of the hot bubble is the densest and least buoyant. A strong buoyancy gradient 
develops which may allow the interior to rise relative to the periphery. In our calculations of 
the mass and area of the hot bubble, we include only the buoyant region of the hot bubble. 
For computational simplicity, we assume that the buoyant region remains intact. Thus, we 
calculate the force on this region and divide by its mass to find its average acceleration. We 
then iterate using small time steps to find the velocity and vertical displacement. 

As the bubble rises, it passes through successively more rarefied ambient gas. In some 
cases, a bubble will rise from its birth neighborhood into the overlaying neighborhood. At 
that point we account for the new environment's effects on the remnant by switching to the 
simulation results for a remnant evolving in the new environment. Thus, we use the first 
part of the simulation for a remnant in the birth neighborhood and the second part of the 
simulation for a remnant in the overlaying neighborhood. We must deal with the possibility 
that the two remnants may have differing lifetimes and hence differing maturation rates. 
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Although the aging process may not be strictly linear, we define a remnant's "maturity" as 
the ratio of the present age to the total lifetime, where the remnant's lifetime is the time 
required for all of the material in a stationary remnant to cool below 2 x 10*^ K (see Paper 
1 for further discussions of remnant death). When we switch from the remnant evolving in 
the lower neighborhood to the remnant evolving in the overlaying neighborhood, we match 
maturities rather than ages. Later in this article, we will perform time integrations. For the 
integrations, we will not count the time gap between the age of the first remnant and the 
age of the second. If we were to count that time, we would overestimate the effects of the 
remnants. Nonetheless, when we allow the remnants to be buoyant, they live longer than 
stationary remnants, as shown in Figure 1. 
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Fig. 1. — In each figure, the dashed hnes track the non-buoyant SNRs and the sohd hues 
track the buoyant SNRs as they rise. In the left and right figures, respectively, the drag 
coefficient, C, which is used to predict the trajectories of the buoyant SNRs, is set to 1.0 and 
0.0, respectively. The height refers to the height of the center of the remnant with respect 
to the Galactic plane. The height gains for this choice of simulation parameters (explosion 
energy = 1.0 x 10^^ ergs, effective magnetic field strength = 5.0 fiG) lie between those for 
the other two sets of simulation parameters used in this paper and described in Section 2.5. 
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2.3. The Ambient Medium in the Thick Disk and Lower Halo 

Both the hydrocode simulations and the buoyancy calculations require estimates of the 
ambient density of material. From Ferriere (1998), the vertical distributions of molecular, 



neutral, warm ionized, and hot hydrogen nuclei are: 

n^{z) = 0.58 e-^^/^^P'^)' cm'^ (4) 

nn{z) = (0.395 e-(^/^27pc)2 ^ q^^q^ g-(./318pc)2 ^ Q ^g^ g-kl/403pc^) ^^-3 (5^) 

n^.(z) = (0.0237 e-l^l/^ + 0.0013 e'l^l/^^" p") cm'^ (6) 

nh{z) = 4.8 X 10-^ e-l^l/^-^'^P'^ cm-^. (7) 



The ratio of helium to hydrogen atoms is approximately 1 to 10, so the total density of atoms 
is assumed to be approximately 1.1 times the density of hydrogen atoms. From Equations 4 
through 7 we see that the H I scaleheight is approximately 130 pc, half of the molecular, 
atomic, and ionized material lies below 2; ~ 90 pc, and 80% of the material lies below the 
scaleheight for isolated supernovae (~ 300 pc, presented below). 

2.4. The Distribution of Isolated Supernova Progenitors 

From Ferriere (1998), the volumetric rate of Type la supernova explosions near the 
Sun's Galactocentric radius is: 

Rateia(^) = 4.0 x lo-^e-l^l/^sspc kpc-V~^- (8) 

The expected rate of isolated Type lb plus Type II supernova explosions is bracketed by 
estimates taken from Ferriere (1998) (lower range) and calculated from McKee & Williams 

(1997) (upper range): 

Rateib+ii(^) = 1.4 to 2.8 x lO-Sg-l^l/^espc kpc^V"^- (9) 

Note that isolated explosions account for only 40% of the Typelb and Type II rate. The 
remaining explosions occur in clusters near the Galactic midplane. They create superbubbles, 
which are well analyzed in Ferriere (1998), but outside our frame of interest. For the following 
calculations, we will use the total isolated SN explosion rate at the Sun's Galactocentric 
radius, in which the isolated Type lb plus Type II rate is taken as the average of the Ferriere 

(1998) and McKee & Wilhams (1997) estimates: 

RatesN(-s) = (4.0 x io-6e-l^l/325pc + 2.I x iQ-Se-l^l/'^^P") kpc^V"'- (10) 

In more conceptual terms, every million years 2 to 3 isolated progenitor stars explode within 
an open-ended column beginning at 2; = 130 pc, extending to 2; = 00, having a cross sectional 
area of 1 kpc^, and residing at the Sun's Galactocentric radius. 



2.5. The SNR Simulation Parameters 



A remnant born at the SN progenitor scale height will be surrounded by a medium 
of density, Uo, equal to 0.1 atoms cm~^. Remnants born at half the scale height will be 
surrounded by ~ 0.3 atoms cm^"^ gas while remnants born at twice the scale height will 
encounter ~ 0.03 atoms cm~^ gas. In order to account for the variation in ambient density 
above the Galactic disk, we performed simulations for the following sample of seven ambient 
densities: 0.5 atoms cm~^ (corresponding to a height of zi — 76 pc), 0.2 (corresponding to a 
height of Z2 = 190 pc), 0.1 (corresponding to a height of = 310 pc), 0.05 (corresponding to 
a height of Z4 = 480 pc), 0.02 (corresponding to a height of ^5 = 850 pc), 0.01 (corresponding 
to a height of Zq = 1300 pc), and 0.005 atoms cm~^ (corresponding to a height of Zj = 
1800 pc). Z7 is roughly 6 times the scaleheight for isolated supernova explosions. The 
n = 0.5 atoms cm~^ simulations were performed for comparison with the others; they will not 
be included in the extra-planar population. Due to uncertainties in the values of the typical 
explosion energy, we performed simulations for both Eo — 0.5 x 10^^ ergs and Eg — 1.0 x 10^^ 
ergs which represent the typical range of estimated explosion energies. At any given height, 
we have approximated the SNRs' environment as being homogeneous with a temperature, 
T, of 1.0 X 10^ K. At this temperature, the hydrogen in this gas is ionized, as is typical of 
the Reynolds layer and of the plasma surrounding and pre-ionized by supernova remnants. 
The turbulent pressure is thought to be small and is not explicitly included in the hydrocode 
simulations. The magnetic and cosmic ray pressures are swept into a single term called the 
nonthermal pressure, Pnt- Observational and theoretical estimates of their values in the thick 
disk are sparse, with the values chosen for the simulations, Pnt — 1800 and 7200 K cm""^, 
being well within the probable range. If B^ff is the effective magnetic field strength, where 
Pnt — -Bgjj/(87r), then the chosen nonthermal pressures correlate to Bgff — 2.5 and 5.0//G. 
We set Beff = 2.5 /jG only when also setting Eg — 0.5 x 10^^ ergs. Thus for every choice 
of ambient density in our sample, we have performed three simulations, the first having 
Eo = 0.5 X 10^^ ergs and -Be// = 2.5;LiG, the second having Eo = 0.5 x 10^^ ergs and 
Beff = 5.0/uG, and the third having Eo = 1.0 x 10^^ ergs and -Be// = 5.0/iG. In all, we report 
on 21 simulations. Their parameters are listed in Table 1. 
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Table 1. Parameters for SNR Simulations 



z 


no 




Beff 


Name 


(pc) 


(atoms cm~^) 


(10^^ ergs) 






76 


0.5 


0.5 


2.5 


halo 65 


;? 


;? 




5.0 


halo 62 






1.0 




halo 64 


190 


0.2 


0.5 


2.5 


halo 59 








5.0 


halo 57 






1.0 




halo 58 


310 


0.1 


0.5 


2.5 


halo 92 








5.0 


halo 91 






1.0 


11 


halo 90 


480 


0.05 


0.5 


2.5 


halo 110 


?) 


11 


)) 


5.0 


halo 111 


)) 


11 


1.0 




halo 109 


850 


0.02 


0.5 


2.5 


halo 100 








5.0 


halo 101 


)) 


11 


1.0 


11 


halo 99 


1300 


0.01 


0.5 


2.5 


halo 103 


)) 


11 


)) 


5.0 


halo 104 




11 


1.0 


11 


halo 102 


1800 


0.005 


0.5 


2.5 


halo 106 


5) 


5) 


11 


5.0 


halo 107 


5) 


5) 


1.0 


11 


halo 105 



- 12 - 



3. Simulation Results 

3.1. Individual Isolated Remnants 

Each computer simulation yielded detailed, comprehensive predictions of the remnant's 
physical state as a function of time. The temperature, pressure, velocity, high ion content, 
high ion absorption line profile. X-ray spectrum, and X-ray luminosity of a sample remnant 
as a function of time are presented in Papers 1 and II. In the present paper, we are focusing 
on the bulk physical and observable effects of the population of isolated extraplanar SNRs 
and so we will extract a complementary set of information from each simulation. 

First, we turn our attention to the cross sectional area covered by supernova remnants. 
When we look up into the Earth's atmosphere we see that some fraction of the sky is covered 
by the population of clouds. Similarly, we can image that some fraction of the Galactic high 
latitude sky is covered by the population of supernova remnants. Determining what this 
fraction is will help us to interpret high latitude observations of hot gas. Theoretically, the 
fraction of sky covered by SNRs above a height Z is RatesN(-2) x (J TTR'^{z,t) dt) dz, 
where RatesN is given in Equation 10 and R{z,t) is the radius of a remnant of age t and 
evolving at height z. In this subsection, we determine the second factor, J 7rR'^{z,t) dt, the 
time integrated area of a SNR evolving at height z. (We will use the integral to calculate 
the fractional area coverage in the following subsection.) We approximate the integral as 
'^7iR'^{z,t) X At. We use up to 26 time slices for each summation and, unless stated 
otherwise, we sum through the end of the remnant's life. The radius is determined from the 
location of the shock front when the remnant is young and from Paper I's bubble boundary 
criterion when the remnant is old enough to have a cool shell. According to this criterion, 
previously shocked material that has cooled to a temperature less than 2 x 10^ K is considered 
to be part of the shell or environmental ISM while hotter material is considered to be part 
of the hot bubble. The time integrated area of each simulated SNR is presented in Table 2. 
We added a column listing the integrated area for the remnants' early stage, the pre-shell 
formation (PSF) phase, which ends when the gas immediately behind the shockfront is no 
longer heated to more than 3 x 10"^ K. This column is identified by the "psf" notation in 
the summation. The pre-shell formation phase is worth attention because during this stage 
the dense gas behind the shockfront is very hot and hence luminous in soft X-rays and UV 
photons. As a result, the remnant could probably be observed and identified as a distinct 
entity during this stage. At the end of this stage, the gas behind the shockfront undergoes 
rapid radiative cooling and so evolves into a cool shell which separates the ever-expanding 
shockfront from the hot, rarefied interior. Although the remaining hot interior (the 'hot 
bubble') can be sufficiently hot and ionized to produce soft X-rays and UV photons, its 
density and therefore its intensity will be less than previously. As a result, an old remnant 
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might not be observed or identified as a distinct entity. Nonetheless, as we will show later, 
the vast majority of the ions (in a time integrated calculation) result from the post-shell 
formation phase. Note that "halo 106" has a larger pre-shell formation area integral than 
"halo 107" and "halo 108" because its shell formed later. 

Let us turn our attention to the volume occupied by SNRs. At any given moment, some 
regions of the Galaxy are occupied by SNRs while others are not. Ignoring overlap, the 
average fraction of space at a height of z which is occupied by SNRs (the "filling fraction" ) 
can be estimated from: RatesN(-2) x (/ ^'n'R^{z,t)dt). In this subsection, we determine the 
second factor in the overall integration, f ^nR^{z,t) dt, the time integrated volume. We 
approximate this integral (which is called the "four- volume" in Slavin & Cox (1993)) as 
^^7rR^{z,t) At. Table 2 presents the time integrated volumes of individual simulated 
remnants. These values will be used in the following subsection to calculate the filling 
fraction of SNR gas. 

In preparation for calculating the average 0+^, N+^, and C"*"^ column densities in the 
following subsection, we calculate the time integrated number of ions per SNR here. The 
integrals are J nion{f, t) 47rr^ dr dt, where nion{f, t) is the volume density of the chosen ion 
at the radius r from the remnant's center and at a time t. The integrals are approximated by 
summations with respect to time and radius and presented in Table 3. Note that C"*"^, and 
to a lesser extent, N"*"^ can be produced via photoionization or collisional ionization. Here, 
we report on only the collisionally ionized C+'^ and N"*"^. As in Table 2, additional columns 
are provided for the integrals through the pre-shell-formation ("psf") stage of evolution. 

Similarly, we calculate the energy released in O"*"^, N"*"^, and C^^ resonance line doublet 
photons (0+^: 1032, 1038 A; N+^: 1239, 1243 A; C+-^: 1548, 1551 A). The energy released 
by a single SNR in a particular doublet is approximated as the doublet's luminosity summed 
over time, At. See Table 4. We also calculate the energy released in 1/4 keV soft 

X-ray photons. The time integrated luminosity is approximated as keV At. The soft 

X-ray photons are scattered across a spectrum. So, in order to compare our predictions 
with observations and other calculations, we convolve the spectra with the i? 05^ T response 
functions for the Rl and R2 bands and report the results in units of ROSATRl + R2 counts 
cm~^. See Table 5. As in the other tables, a column is provided for the contribution from the 
pre-shell formation phase, when the remnants are young, bright, and most easily identified 
as individual objects. 
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Table 2. 



Simulated 
Remnant 


JJo'^ 7ri?2 X At 
(kpc^ yr) 


J2 T^R^ X At 

(kpc^ yr) 


-7ri?3 X At 
(kpc'^ yr) 


\t^R^ X At 
(kpc^ yr) 


Uo — 0.5 cm~^ 
halo 65 
halo 62 
halo 64 


16 
35 
72 


25,000 
18,000 
29,000 


0.35 
0.77 
2.2 


1500 

980 

1900 


Uo = 0.2 cm~^ 
halo 59 
halo 57 
halo 58 


81 
81 
110 


61,000 
31,000 
58,000 


2.7 
2.7 
4.2 


4700 
2100 
4800 


no = 0.1 cm~^ 
halo 92 
halo 91 
halo 90 


360 
220 
300 


110,000 

49,000 
91,000 


18 
10 
15 


10,000 

3800 
8700 


rio = 0.05 cm~^ 










halo 110 
halo 111 
halo 109 


960 
500 
1300 


170,000 
75,000 
140,000 


63 
29 
100 


19,000 
6600 
15,000 


no = 0.02 cm-''^ 










halo 100 
halo 101 
halo 99 


740 
750 
1000 


300,000 
140,000 
260,000 


53 
54 
83 


40,000 
15,000 
34,000 


Uo = 0.01 cm"^ 










halo 103 
halo 104 
halo 102 


1000 
1000 
1400 


480,000 
250,000 
460,000 


84 
87 
130 


72,000 
31,000 
69,000 


Uo = 0.005 cm-3 










halo 106 
halo 107 
halo 105 


6800 
1400 
1800 


800,000 
500,000 
900,000 


940 
140 
210 


140,000 

73,000 

160,000 
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Table 3. 



Simulated 
Remnant 


Y^psf 0+5 At 

(0+5 ion yr) 


E 0+5 At 
(0+5 ion yr) 


(N+'* ion yr) 


J2 N+* At 
(N+4 ion yr) 


(C+^ ion yr) 


X; c+3 At 

(C+^ ion yr) 


halo 65 
halo 62 

Viol/-, 


2.7 xl058 
2.1 xl059 

O.U XiU 


9.2 xlO^o 
8.5 x 10*50 

i.y X lu 


1.8 xl057 
1.6 xl058 

4.0 XiU 


1.0 xlO^o 
9.0 xl059 
i.y XiU 


1.4 xl057 
6.6 xl058 

O.i XiU 


2.4 xlO^o 
1.6 X 10*50 

o.z XiU 


halo 59 
halo 57 
halo 58 


2.7 xl059 
3.2 xl059 
2.0 xl059 


1.9 x 10*51 
2.0 x 10*51 
4.2 x 10*51 


1.1 Xl058 

1.2 xl05» 

1.2 Xl058 


2.2 X 10*50 
1.6 X 10*50 

3.3 X 10*50 


8.6 xl057 
1.5 xl058 
1.5 xl058 


4.6 X 10*50 

2.7 X 10*50 
5.4 X 10*50 


halo 92 
halo 91 
halo 90 


1.6 xlO'^o 
6.4 xl059 

8.1 Xl059 


3.4 x 10*51 
3.6 X 10*51 

7.5 X 10*51 


9.9 xl058 
6.6 xl058 
3.4 xl058 


3.8 X 10*50 
2.6 X 10*50 
5.0 X 10*50 


1.6 Xl059 

6.1 xl058 

3.7 xl058 


7.5 X 10*50 
4.7 X 10*50 
9.0 xlO^o 


halo 110 
halo 111 
halo 109 


2.8 xlO^" 

1.4 xlO*^" 

4.5 xlO^o 


6.4 X 10*51 
6.0 X 10*51 
1.3 X 10*52 


2.9 xl059 
8.7 xl058 
3.0 xl059 


6.7 X 10*50 
4.3 X 10*50 
8.3 X 10*50 


4.7 xl059 
2.4 xl059 
5.1 xl059 


1.3 xlO*5i 
8.9 x 10*50 
1.6 X 10*51 


halo 100 

halo 101 
halo 99 


1.3 xlO^° 

1.7 xlO^o 

1.8 xlO^" 


1.4 X 10*52 

1.2 X 10*52 
2.4 X 10*52 


7.7 xl05^ 

1.5 Xl059 
1.1 Xl059 


1.4 X 10*51 

8.9 X 10*50 
1.7 xl0*5i 


9.8 xl05S 

1.8 Xl059 

1.3 xl05" 


2.8 X 10*51 

2.1 xl0*5i 

3.9 X 10*51 


halo 103 
halo 104 

halo 102 


1.3 xlO^o 
1.7 X 10*^0 
1.9 X 10*50 


2.4 X 10*52 
1.9 X 10*52 

3.9 X 10*52 


8.1 xl058 

1.3 Xl059 
1.3 Xl059 


2.3 X 10*51 
1.5 X 10*51 

2.9 X 10*51 


1.2 Xl059 

1.7 Xl059 

1.8 Xl059 


5.2 X 10*51 
3.7 X 10*51 
7.2 X 10*51 


halo 106 
halo 107 
halo 105 


8.3 X 10*50 
2.1 xlO^o 
2.6 xlO^o 


3.9 X 10*52 
3.7 X 10*52 
7.6 X 10*52 


6.7 xl059 

1.5 Xl059 
1.9 Xl059 


3.7 X 10*51 

3.3 xlO^i 

6.4 X 10*51 


2.8 X 10*50 
2.4 xl059 
3.0 xl059 


9.4 X 10*51 

8.6 xlO^i 

1.7 X 10*52 
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Table 4. 



Simulated 
Remnant 


Ya)'''^ ^o+^ Af 
(ergs) 


(ergs) 


Eo"'^ iw+i Ai 
(ergs) 


E At 
(ergs) 


Eo*'^^c+3 At 
(ergs) 


At 

(ergs) 


halo 65 
halo 62 

Viol/-, 


5.6 xl0^7 

3.2 xlO*^ 

1 /I V 1 n49 


1.0 xlO*9 
1.0 xlO*9 

o.y xiu 


5.4 xl0^6 
2.4 xlO^^ 


1.8 xlO*^ 

7.4 xlO''^ 

7 1 V 1 n48 
( . i X iU 


4.6 xl0^6 
2.4 xlO^^ 

o.U XlU 


7.4 xlO** 
7.4 xlO*^ 

1.4 XiU 


halo 59 
halo 57 
halo 58 


3.3 xio-^s 
3.9 xlO*^ 
2.1 xlO^s 


1.4 xW^ 
1.6 xlO*9 
2.8 xlO'^s 


1.6 xlO^'' 

1.7 xlO^^ 
1.4 xlO''^ 


1.8 xlO^s 
1.6 xlO** 
3.4 xlO'^s 


1.3 xl047 
2.2 xlO^^ 
1.8 xlO''^ 


5.0 xW* 
4.4 xlO** 

1.1 xlO^s 


halo 92 
halo 91 
halo 90 


9.4 xlO"** 
3.4 xlO^s 

4.6 xlO*^ 


1.6 xlO^^ 
1.5 xW^ 
3.0 xl0*9 


7.6 xlO^^ 

5.2 xlO^^ 

2.3 xlO^^ 


1.9 xlO^s 
2.2 xlO^s 

3.2 xlO** 


1.1 xlO^s 
5.8 xlO-*^ 
2.5 xlO^^ 


4.6 xlO^* 
4.2 xlO^s 

9.7 xlO** 


halo 110 
halo 111 
halo 109 


7.6 xlO"*** 
3.6 xlO^s 
1.2 xlO'^s 


1.4 xlO^^ 

1.5 xlO''^ 
3.4 xlO^^ 


1.1 xlO^s 
2.8 xl047 
1.0 xlO^^ 


2.2 xlO^* 
1.6 xlO^* 
3.5 xW* 


1.6 xlO''^ 
9.0 xlO*^ 
1.4 xlO^^ 


4.2 xlO^* 
4.2 xlO^* 
6.4 xlO^* 


halo 100 

halo 101 
halo 99 


1.2 xlO^^ 

1.4 xlO'** 
1.4 xlO''* 


1.2 xlO*9 

1.9 xlO'*'^ 
3.5 xl049 


9.1 xlO^^ 
1.7 xl047 
1.0 xl047 


1.6 xlO*^ 
2.1 xlO^s 
3.5 xlO^s 


1.3 xlO^^ 
2.6 xlO''^ 

1.4 xlO"^ 


2.7 xlO*^ 

3.3 xlO^* 
5.0 xlO'*^ 


halo 103 
halo 104 

halo 102 


4.0 xlO''^ 

5.1 xlO*^ 
5.1 xlO^^ 


1.2 xl049 
2.0 xlO*9 

3.8 xlO^^ 


3.1 X 10^*5 

5.0 xl0^6 

4.1 xlO"*^ 


1.6 xlO''* 
1.9 xlO** 

3.6 xlO^s 


5.7 xlO^'s 
8.6 x 10^*5 
7.1 xlO'"' 


3.0 xlO''* 
3.0 xlO** 
5.6 xlO^s 


halo 106 
halo 107 
halo 105 


1.1 xlO^'s 

2.2 xlO*7 
2.4 xlO*^ 


1.5 xlO'^s 
2.2 xlO*9 
4.5 xlO*9 


1.3 xlO''^ 
2.1 xl0^6 
2.1 xl0^6 


1.9 xlO'^s 

2.3 xlO** 

4.4 xlO*^ 


6.3 xlO^^ 

4.1 xl0^6 

4.2 xl0^6 


4.6 xlO'^s 
4.3 xlO** 
8.8 xlO** 
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Table 5. 




Simulated 
Remnant 


Ylo^'^ -^1/4 keV X '^t 

(ROSAT Rl + R2 counts cm^) 


keV X 

(ROSAT Rl + R2 counts cm^) 


halo 65 
nalo uz 
halo 64 


6.9 xlO^'^ 
2.2 xlO''" 
9.9 xlO^^ 




2.4 xlO^o 

0. i X iU 

3.1 xlO^o 


halo 59 
halo 57 
halo 58 


8.4 xlO^^ 
8.0 xlO^^ 
1.4 xlO^o 




1.7 xlO^o 

i. ( X iU 

4.0 xlO<^o 


halo 92 
naio yi 
halo 90 


9.5 xlO^^ 

7.2 xlO^'^ 

1.3 xlO^o 




1.2 xlO^o 

i.o XiU 

3.0 xlO^o 


halo 110 

1, „ 1^ 111 
nalo iii 

halo 109 


7.0 xlO^'^ 

5.1 xlO^^ 
1.3 xlO^" 




9.0 xlO^^ 

1.1 XiU 

2.5 xlO^o 


halo 100 
nalo iui 
halo 99 


2.6 xlO^^ 
2.3 xlO^^ 
4.1 xlO^^ 




7.2 xl0^9 

L.Z XiU 

2.9 xlO^o 


halo 103 
naio iU4 
halo 102 


1.5 xlO^^ 
1.3 xlO^'^ 
2.3 xlO^^^ 




5.9 xlO^^ 

i.o XiU 

3.5 xlO^° 


halo 106 
halo 107 
halo 105 


2.1 xlO^'^ 
7.4 xlO^s 

1.2 xl0^9 




5.8 xlO^^ 
1.5 xlO^o 
3.7 xlO^o 
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3.2. SNR Population: Area Coverage and Volume Occupation 

In this and the following two subsections, we predict the effects of the population of 
isolated extra-planar supernova remnants. For example, we will calculate the fraction of the 
high latitude sky's area that is covered by supernova remnants residing above a height, Z, 
of 130 pc. The necessary integral, presented in Section 3.1, will be approximated by the 
summation Y1'z(^^'^^S'n{z) x '^7rR'^{z,t)At). The first step in handhng the gradients in 
ambient density and SN rate with respect to z is to segment the thick disk and halo into six 
plane parallel slabs whose boundaries reside at the midpoints between zi, Z2, Z3,Z4, z^, zq, 
and zr. The top slab extends to z = oo. (This parceling scheme begins sX z = 130 pc, the 
scaleheight for the Galactic H I layer, thus excludes the ^i, no = 0.5 cm~^ simulations which 
had been created for comparison with the other simulations.) 

We then integrate RatesN(-2), the isolated SN rate per unit volume and time given in 
Equation 10, from the lower to the upper boundary of each slab. We multiply each slab's 
rate integral by the time integrated cross sectional area for a single SNR residing in that slab 
tabulated in Table 2. This act implicitly treats all SNRs residing within a slab as if they 
are residing at the representative height of the slab (i.e. Z2i z^, etc.). For each of the slabs. 
Table 6 lists the integrated SN rates and the fraction of area covered by the population of 
SNRs if the remnants are stationary. Next, we calculate the effects of the entire population of 
isolated SNRs born above 130 pc. We sum the fractional area coverages of the six slabs. For 
stationary SNRs, we use the area coverage values listed in Table 6. For buoyant SNRs, we 
adjust the time integrated cross sectional areas to account for the time spent in the various 
slabs. The results for both buoyant and non-buoyant SNRs are tabulated in Table 7. Note 
that we maintain accuracy to several significant digits during all calculations, but round to 
two significant digits when reporting results. 

Table 6 implies that about 1% of the high latitude sky is covered by young, presumably 
bright and recognizable SNRs (the pre-shell formation remnants). Depending on the assumed 
simulation parameters and drag coefficient, 30% to 90% of the sky is covered by the hot gas 
in SNRs of any age^. Although the older SNRs that account for most of the area coverage 
are not as recognizable as their younger cousins, they still harbor large stores of 0+^, N"*"^, 
and which are observationally important. 

The "filling fraction" , the fraction of volume at a given height that is filled by extraplanar 
supernova remnants is simply: RatesN(^) x ^ ^nli^{z, t) At. We calculate the filling fraction 



^The cross section of a SNR residing at a given height may overlap with the cross section of a SNR residing 
at another height. Therefore the probabiHty that one or more SNR He along a sightline perpendicular to the 
Galactic plane is somewhat smaller than the sums presented. 
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at the representative height of a slab (i.e. Z2, Z3, z^, Z5, zq, or zj) by calculating the SN 
explosion rate at the chosen height and then multiplying by the time integrated volume 
occupation of a single SNR born at the chosen height (taken from Table 2). The resulting 
filling fractions for non-buoyant SNRs are listed in Table 6. At all examined heights, isolated 
non-buoyant SNR bubbles fill less than 10% of the volume. If the remnants are buoyant, 
then the filling fractions in the first couple of slabs would be even smaller, while the minute 
filling fractions in the upper slabs would be somewhat larger. The smallness of our estimates 
echo those from other studies, such as Ferriere (1998). On account of their small volume 
occupations, it is very unlikely that hot supernova remnant bubbles will collide with each 
other. Interestingly, the volume filling fraction reaches its maximum near the SN scaleheight. 
This is because the SN rate decreases as a function of z, the remnants' size and longevity 
increase as a function of z, and the maximum product occurs around z — 300 pc. 
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Table 6. 



k5lab 


Height 


J oW Kate az 


fab Area 


lotal Area 


Fbr Volume 


lotal Volume 






(kpc yr ) 


Coverage [/o) 


L-overage [/o) 


r iliecl [ /o) 


r illea y/o) 


Eo = 0.5 X 10^1 ergs 














B.ff = 2.5 H G 














Yi = 0.2 a,toin cm ^ 


130 to 250 


1.5 xlO~^ 


0.012 


8.9 


0.0034 


6.0 


71 — 0.1 cbtoin cm ^ 


250 to 390 


1.15 xlO~^ 


0.041 


12 


0.014 


7.9 


Ti 0.05 atom cm ^ 


390 to 670 


1.0 xlO~^ 


0.099 


18 


0.024 


7.2 


72 = 0.02 a.tom cm ^ 


D^U T^O i±UU 


A Q V 1 0""^ 
X iU 


U.UoD 


±0 




A 7 


71 = 0.01 atom cm""^ 


1100 to 1600 


1.2 xlO"'' 


0.012 


5.6 


0.0020 


1.7 


n = 0.005 atom cm~^ 


1600 to oo 


2.5 xlO"* 


0.017 




0.00028 


0.041 


Sum 


1.30 to oc' 


4.3 xlO-*' 


0.22 


61 






Eo = 0.5 X 10''^ ergs 














Be// = 5.0 fi. G 














n = 0.2 atom cm^ "' 


130 to 250 


1.5 xlO"^ 


0.012 


4.6 


0.0034 


2.7 


n = 0.1 atom cm~'^ 


250 to 390 


1.2 xlO-s 


0.026 


5.6 


0.0078 


3.0 


n = 0.05 atom cm~^ 


390 to 670 


1.0 xlO-s 


0.051 


7.7 


0.011 


2.5 


n = 0.02 atom cm~^ 


670 to 1100 


4.8 xlO-'^ 


0.036 


6.8 


0.0063 


1.7 


n = 0.01 atom cm~^ 


1100 to 1600 


1.2 xlO-''' 


0.012 


3.0 


0.0021 


0.73 


n = 0.005 atom cm~* 


1600 to oo 


2.5 xlO-8 


0.0035 


L3 


4.2 x 10-5 


0.022 


Sum 


130 to oo 


4.3 xlO-6 


0.14 


29 






Eo = 1.0 X 10^1 ergs 














Be// = 5.0 Ai G 














71 = 0.2 atom cm^'* 


130 to 250 


1.5 xlO-s 


0.016 


8.4 


0.0053 


6.0 


n = 0.1 atom cm^'' 


250 to 390 


1.2 xlO-i^ 


0.034 


11 


0.012 


6.8 


n = 0.05 atom cm~^ 


390 to 670 


1.0 xlO"'^ 


0.14 


14 


0.038 


5.8 


n = 0.02 atom cm~^ 


670 to 1100 


4.8 xlO-'' 


0.048 


12.5 


0.0097 


3.9 


n = 0.01 atom cm~^ 


1100 to 1600 


1.2 xlO-'^ 


0.016 


5.4 


0.0031 


1.6 


n = 0.005 atom cm~^ 


1600 to oo 


2.5 xlO-8 


0.0046 




6.3 X 10-5 


0.049 


Sum 


21 to oo 


6.4 xlO-8 


0.25 


54 
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Table 7. Area covered by non-buoyant or buoyant SNR bubbles. Case 1 assumes 
Eg = 0.5 X 10^-*^ ergs and B^ff = 2.5 jj,G. Case 2 assumes Eo = 0.5 x 10^^ ergs and 



Beff = 


5.0 /iG. Case 3 


assumes = 1.0 x 10^^ 


ergs and Beff = 5.0 /iG 


Case 


Coverage (%) 


Coverage (%) 


Coverage (%) 




(stationary) 


(drag coefficient = 1) 


(drag coefficient = 0) 


1 


61 


65 


87 


2 


29 


29 


32 


3 


54 


55 


67 
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3.3. The C+3, N+^, and Predictions for the Population of SNRs 

Here, we estimate the time and space averaged column densities of collisionally ion- 
ized O^^, N"*"^, and C^^ ions along high latitude lines of sight due to the population of 
isolated extraplanar SNRs. The column density of an ion residing above a height of Z is 
J^(RatesN(^) X / ^mrSj^t) 47rr^ dr dt) dz. We replace the integrals with summations. 
For each slab, we integrate the supernova explosion rate with respect to z from the bottom 
to the top of the slab (see Table 6), and multiply by the 0+^, N+^, and C"*"^ contributions 
made by a SNR born within the slab (see Table 3 for stationary SNRs). When calculating 
the O"*"^, N"*"^, and contributions made by buoyant SNRs, we take into account the ex- 
tended lifetime of the SNR and the time spent in more rarefied environments. We then sum 
across all slabs above z = 130 pc. Table 8 tabulates the resulting average column densities, 
for remnants that are assumed to be buoyant as well as those that are not. For example, the 
high latitude sky average 0+^ column density ranges from 2.3 to 6.3 x 10^^ cm~^, depending 
on the simulation parameters and drag coefficient. This range may appear to be quite com- 
pact. Moderately buoyant SNRs would produce only 3 to 20% more O^^ than stationary 
SNRs, while unfettered SNRs would produce 35 to 110% more O"*"^ than stationary SNRs. 
For identical drag coefficients. Cases 1, 2, and 3 vary by a factor of 2.1 to 2.4, primarily 
due to the factor of 2 variation in assumed explosion energy. Interestingly, increasing the 
nonthermal pressure by a factor of 4 (from Case 1 to Case 2) slightly decreases the net 
O"*"^ contribution. The reason for this phenomenon is that remnants experiencing greater 
nonthermal pressure are smaller, therefore denser, therefore more luminous, and therefore 
shorter lived. Note that Table 8's values were found by averaging over many SNRs and even 
parts of the sky lacking SNRs. For predictions of 0+^, N+^, and C+^ column densities within 
individual SNRs, as a function of age and radius, see Paper I. 

The height distribution of the 0+^-rich gas in extraplanar SNRs is illustrated in Figure 2. 
For all choices of simulation parameters and drag coefficients, the z — 390 to 670 pc slab 
makes the largest contribution to the sky averaged 0+^ column density. In Section 4.2, we 
will add the local 0+^ distribution and compare the result with Copernicus, ORFEUS, and 
FUSE observations. 

Figure 3 displays the distribution of buoyant velocities for the population of supernova 
remnants. In making these plots, we assumed that the drag coeffient is or 1. No plot is 
needed for the stationary SNRs, because their buoyant velocities are km s~^. The plotted 

buoyant velocities are directed away from the Galactic plane. The first step in calculating 
the distributions was to calculate each remnant's buoyant velocity as a function of time 
by integrating the buoyant acceleration with respect to time using small time steps. We 
assumed that buoyancy does not tear the remnants apart, so the calculated buoyant velocity 
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was assumed to pertain to the entire supernova remnant. Remnants born at different heights 
above the plane have different buoyancy histories. Thus, the final distributions of velocities 
were found by making weighted compilations of the remnants' buoyancy histories, where 
the weightings were proportional to the supernova explosion rate in the given slab. Case 1 
{Eo = 1.0 X 10® ergs, -Be// = 5fiG) remnants are the most buoyant, probably because the 
small assumed ambient nonthermal pressure allows the hot bubbles to expand to greater 
volumes, lowering the internal density and raising their buoyancy. Case 2 remnants are the 
least buoyant for the converse reason. With drag coefficients of and 1, we only see the 
0"''^-rich gas moving away from the plane; we do not see it stall and begin to fall back 
toward the plane. However, we saw the 0"'"^-rich gas rise and fall in preliminary estimates 
using larger drag coefficients. If there is moderate drag (C = 1) then the typical buoyant 
velocity will hover around 20 km sec~^. If there is no drag (C = 0), then the typical buoyant 
velocity increases to ~ 70 to ~ 100 km sec~^. If we were to observe such remnants, relative 
motions within the Galaxy would probably increase the spread of observed velocities, as 
would motion within the remnants themselves (see Paper I for estimates). Observations 
pointed at non-zenith directions record only the radial velocity, a trigonometric fraction of 
the velocity in the vertical direction. 
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Fig. 2. — The O^^ column density expected to reside on vertical sightlines extending from 
the Earth to given heights above the plane, as a function of height above the plane. Various 
choices of buoyancy are marked by color. Green denotes a population of stationary SNRs. 
Blue denotes a population of moderately buoyant SNRs (drag coefficient = 1). Purple 
denotes a population of unrestrainedly buoyant SNRs (drag coefficient = 0). Independent of 
the choice of drag coefficient, the z = 390 to 670 pc slab makes the greatest contribution to 
the total column density through the halo. Note that the 0+^ appears to begin at 2; = 190 pc, 
only because we set off the lowest SN explosions at the representative height for the lowest 
included slab, which is Z2 = 190 pc. As in Figure 1, the simulated SNRs have an explosion 
energy of 1.0 x 10^^ ergs and ambient effective magnetic field of 5.0/xG (i.e. Case 3). 
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Fig. 3. — Distribution in buoyant velocity among the population of extraplanar supernova 
remnants, assuming moderate drag (C = 1, left figure) or no drag (C = 0, right figure). In 
addition to buoyant bulk velocity directed away from the Galactic plane, the young remnants 
will also exhibit expansion velocities. Detailed predictions of the expansion are provided in 
Paper I. 
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The proceedure for calculating the high ion resonance line intensities is very similiar to 
that for calculating the column densities. Theoretically, the average intensities are given by 
J^(RatesN(2;) x f L{z,t)dt) (i2;/(47rsteradians). We replace the integrals with summations. 
The resulting intensities are listed in Table 9. Allowing unrestrained buoyancy increases the 
sky-averaged intensities by only 20% or less, thus much less significantly than it affects the 
high-stage ion column densities. Buoyancy mostly affects the very late-term evolution of 
the remnants by extending their lives and elevating them into more rarefied surroundings. 
However, old remnants are dim and moving them to more rarefied surroundings makes them 
even dimmer. Therefore, these two influences on net photon production counteract each 
other. For the same reason, we see that buoyancy has almost no effect in the following 
calculations of 1/4 keV X-ray intensity. 
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Table 8. The model population of SNRs in the thick disk produce the following 
sky-averaged O^"'"', N^'^, and C+'^ column densities. Case 1 assumes Eo = 0.5 x 10^^ ergs 
and Bgff = 2.5 /xG. Case 2 assumes Eo = 0.5 x 10^^ ergs and B^ff = 5.0 /iG. Case 3 
assumes Eo = l.O x 10^^ ergs and Bgff = 5.0 /iG 



Ion and Case 


Column Density 


Column Density 


Column Density 




(ions cm~^) 


(ions cm~^) 


(ions cm~^) 




(stationary SNR) 


(drag coefficient = 1) 


(drag coefficient — 0) 


0+5 








1 


2.5 X 10^3 


2.7 X 10^3 


4.0 X 10^3 


2 


2.3 X 10^3 


2.4 X 10^3 


2.6 X 10^3 


3 


4.8 X 10^3 


5.0 X 10^3 


6.3 X 10^3 










1 


2.6 X 10^2 


2.9 X 10^2 


4.2 X 10^2 


2 


1.7 X 10^2 


1.8 X 10^2 


2.1 X 10^2 


3 


3.4 X 10^2 


3.5 X 10^2 


4.7 X 10^2 










1 


5.3 X 10^2 


6.3 X 10^2 


1.1 X 10^3 


2 


3.7 X 10^2 


3.8 X 10^2 


5.0 X 10^2 


3 


7.0 X 10^2 


7.4 X 10^2 


1.2 X 10^^ 
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Table 9. The model population of SNRs in the thick disk produce the following 
sky-averaged O^''"^ N^"*, and C"*"^ resonance line doublet intensities. Case 1 assumes 
Eg = 0.5 X lO^-*^ ergs and B^ff = 2.5 /xG. Case 2 assumes Eo = 0.5 x 10^^ ergs and 
Beff = 5.0 iiG. Case 3 assumes Eo = 1.0 x 10^^ ergs and B^jf = 5.0 /xG 



Ion and Case 


Intensity 


Intensity 


Intensity 




(photon s~^ cm~^ sr~^) 


(photon s~^ cm~^ sr~^) 


(photon s~^ cm~^ sr~^) 




(stationary SNR) 


(drag coefficient = 1) 


(drag coefficient = 0) 


0+5 








1 


840 


850 


920 


2 


940 


960 


1000 


3 


1800 


1900 


2100 










1 


130 


140 


150 


2 


130 


130 


140 


3 


240 


250 


270 










1 


390 


390 


410 


2 


370 


370 


380 


3 


780 


790 


830 
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3.4. The X-ray Predictions for the Population of SNRs 

Following the procedures outlined in the previous two subsections, we estimate the 
average ^ keV surface brightness due to the ensemble of extraplanar supernova remnants 
(see Table 10). The average 1/4 keV countrate is 120 to 320 xlO^^ counts s^^ arcmin^^, 
depending on the choice of explosion energy, magnetic pressure, and to a minute extent, drag 
coefficient. These estimates would increase by about 50% if we were to add the zi SNRs 
between z — 21 and 130 pc. The predictions are insensitive to buoyancy for the same reason 
as the high ion intensities are little affected by buoyancy - the longer lifetimes of buoyant 
remnants are counteracted by their lesser luminosities in old age. 

Because we have averaged the SNR populations' brightness over the entire high latitude 
sky, the resulting average countrate is less than the countrate from individual remnants, 
especially young remnants. Roughly half of the soft X-rays are emitted during the bright, 
but brief stage before the remnant forms a cool shell. The remaining half are emitted 
later, when the old, dim remnants are difficult to identify. Observationally, the collective 
emission from unidentifiable remnants covering half of the sky would be construed as an 
X-ray background, punctuated by a few bright arcs tracing the limbs of young remnants. 
We expand on this description in the following subsection. 
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Table 10. Sky- Averaged 1/4 keV Countrates (in the ROSAT Rl + R2 bands) from the 

Population of SNRs above z — 130 pc. Case 1 assumes Eo — 0.5 x 10^^ ergs and 
Sg// = 2.5 //G. Case 2 assumes — 0.5 x 10^^ ergs and 5e// — 5-0 //G. Case 3 assumes 

= 1.0 X 10^^ ergs and 5^// = 5.0 //G 

Case Countrate Countrate Countrate 

/-j^g— 6 counts \ flO~^ counts \ flO~^ counts \ 

^ s aJTcmin^ ' ^ s arcmin^ ^ s arcmin^ 

(stationary SNR) (drag coefficient = 1) (drag coefficient = 0) 



1 

2 
3 



120 
130 
310 



120 
130 
310 



120 
130 
320 
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3.5. Spatial AppeEirance 

The high latitude sky (outside of superbubbles) should be comprised of three types 
of regions, those having bright young remnants seen in projection (~ 0.2% of sky), those 
having dimmer older remnants seen in projection (~ 30 to ~ 90% of sky), and those having 
no remnants (~ 10 to ~ 70% of sky). The young remnants should appear edge brightened 
in I keV X-ray emission, 0+^, N+^, and emission, and in numbers of 0+^, N+^, and 
ions. For example, the limb of a pre-shell formation remnant evolving in an n = 0.01 cm~^ 
environment emits more than 3000 xlO^^ | keV counts s~^ arcmin"^. It would be easily 
observed. Remnants evolving nearer to the plane (i.e. at 2; < 1300 pc) would be even 
brighter. 

As the edge of its hot bubble cools below ~ 10^'^ K a remnant loses its sharply edge- 
brightened look. If becomes slightly edge brightened and then centrally brightened (see 
Papers I and II for radial profiles of sample SNRs). Its total luminosity plummets. Our 
sample SNR's ^ keV surface brightness drops to < 500 x 10~^ counts s^^ arcmin^^. An 
old remnant should be difficult to recognize, considering that it may be outshone by the 
unobscured Local Bubble (~ 300 to ~ 800 xlO~® counts s"^ arcmin^^) combined with 
the extragalactic background (~ 400 xlO~^ counts s~^ arcmin"^, de-absorbed). If an old 
remnant were to be detected in the X-ray regime, it should exhibit a soft spectrum and strong 
recombination edges (see Paper II for spectra of sample SNRs) . Rather than being stars of 
the show, these remnants would create mottled soft X-ray and high stage ion backgrounds. 
Note that if O^^ column density and 1/4 keV maps were made from the SNR simulations, 
they would probably look somewhat dissimilar because of the significantly different time 
evolution in these observables. 



4. Compcirison between Simulations and Observations 

4.1. Soft X-rays 

Early, modest angular resolution surveys (the Wisconsin Rocket Program, McCammon 
et al. 1983, McCammon et al. (1983), SAS 3, Marshall & Clark (1984) and the HEAO 1 A-2 
Low Energy Detectors, Garmire et al. (1992)) detected bright soft X-ray (~ \ keV) emission 
from every direction on the sky. The maps show considerable structure and identifiable 
features such as Loop I, the Eridanus Soft X-ray Enhancement, and the Monogem Ring. 
The poles are brighter than the equator and the northern polar region is brighter than the 
southern polar region. Most of the observed soft X-rays originate in the Local Bubble (a 
~ 60 pc radius region surrounding the sun) and the Milky Way's halo. Extragalactic sources 
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produce a fraction of the average observed total (McCammon et al. 1976; Snowden & Pietsch 
1995). Small contributions to the flux observed at low latitudes are made by pre-radiative 
phase remnants and Galactic point sources. (Rosner et al. 1981). 

The halo and extragalactic soft X-rays are attenuated by the layer of diffuse neutral 
gas in the disk, so that only ~ 2/3rds reach the solar system. In directions with especially 
opaque clouds, distinct "shadows" appear in the X-ray maps. It was by observing shad- 
owed regions with the Rontgcnsatallit Position Sensitive Proportional Counter (or ROSAT 
PSPC), that researchers first conclusively demonstrated the emissivity of the distant (halo 
plus extragalactic) component (Snowden et al. 1991; Burrows & Mcndcnhall 1991; Kcrp et 
al. 1993; Snowden et al. 1994; Wang & Ye 1996). More recently, Snowden et al. (1998, 
2000); Kuntz & Snowden (2000) extended this type of analysis to the entire sky. Their 
estimates for the surface brightness of the de-absorbed distant component range from ~ 400 
to 3000 X 10~^ ROSAT 1/4 keV counts s"^ arcmin"^, depending on the direction. The high 
end of the quoted range is misleading, because it refers to the flux from anomalous objects, 
such as the Draco region, and intermediate distance objects, such as the North Polar Spur 
and the Eridion Superbubble, which the analysis technique partitioned between the near and 
far components. Due to the North Polar Spur and Draco, the map of the northern sky is 
brighter and more complex than that of the southern sky. The southern sky provides only 
~ 400 to 1000 X 10~^ counts s~^ arcmin"^, (unabsorbed) and has an average count rate for 
b < —65° of 810 X 10~^ counts s"^ arcmin"^. Because the southern halo is far less "polluted" 
by anomalous regions, this paper takes it as the standard. Next, the extragalactic compo- 
nent must be subtracted from the distant component to yield the halo flux. Prom shadowing 
studies of nearby spiral galaxies, the 1/4 keV extragalactic intensity has been estimated to 
be about 30 keV cm^^ s^^ sr^^ kcV~^, equivalent to 400 x 10~^ counts s^^ arcmin^^ in the 
ROSAT \ keV band (see Snowden & Pietsch (1995); Barber et al. (1996); Cui et al. (1996)). 
The extragalactic surface brightness should be fairly smooth and isotropic. Subtracting its 
flux leaves only ~ 400 x 10~^ counts s^^ arcmin"^ (unabsorbed), which is attributable to 
the Milky Way's southern halo. 

In comparison with our simulation results, 310 xlO^® counts s^^ arcmin"^, or 3/4 of 
the observed emission could come from the population of extraplanar SNRs if their explosion 
energies are 1 x 10^^ ergs, their surrounding effective magnetic field is 5 /xG, and they are 
stationary. Our estimate rises to 400 xlO~^ counts s~^ arcmin"^, or 100% of the observa- 
tionally determined surface brightness if we use McKee & WiUiams (1997) 's isolated massive 
progenitor SN rate rather than averaging it with Ferriere (1998) 's rate. Assuming the more 
modest rate and the full range of explored parameters yields estimated brightnesses of 120 
to 320 xlO~^ counts s~^ arcmin"^, or 30% to 80% of the observationally determined fiux. 
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Other important characteristics of the X-ray emitting gas include its vertical extent, 
degree of spatial non- uniformity, and temperature. Determining the height of the X-ray 
emitting plasma is difficult, but shadowing by low, intermediate, and high velocity clouds 
provides some indication. After comparing the 1/4 keV X-ray emission with the H I column 
densities in and around the M complex, Herbstmeier et al. (1995) concluding that nearly 
all of the distant emission originates beyond the low velocity clouds {zdoud ~ 200 pc), 
while some originates below and some originates beyond the high velocity clouds (1.5 kpc ~ 
^cioud ~ 4.4 kpc). Hence the distribution of extraplanar SNRs is in rough agreement with 
the distribution of X-ray emitting gas. 

The eye can easily identify patcliincss in the 1/4 keV survey maps of Snowden et al. 
(1998). Quantitatively, the halo emission varies by a factor of 2 over angular scales of ~ 20° 
(Snowden et al. 1997). On a smaller scale. Burrows & Mendenhall (1991) found that the 
emission behind the Draco cloud varies by 1/3 from the northern to the southern half of the 
region (~ 1°). In the Milky Way's south, the flux is arranged in a fringe of bright regions 
northward of 6 ~ —45° (which may be artifacts of the deabsorption analysis), a mottled 
background, and a couple of noticeably bright regions below / ~ —45° (i.e. the regions 
around / = 247°, b = —64° and / = 215°, b = —68°). We think that these two bright regions 
must be emitting features rather than artifacts of the analysis technique because there is 
no significant spatial correlation between the distant emission map and the absorption map 
used in Snowden et al. (1997) 's analysis. The shapes of these bright features can be better 
seen in the less smoothed total observed surface brightness map of Snowden et al. (1997). 
This map shows that the larger feature, the one around I — 247°, b — —64° has a fairly crisp 
outer edge and a crescent, edge brightened shape. 

Approximately 1 bright, pre-shell formation phase SNRs is expected to reside in the 
halo of each hemisphere above &~|50°| and z ~ 130 pc. The feature around I — 247°, 
b — —64° is a reasonable candidate for a pre-shell formation phase halo SNR. The Dickey & 
Lockman (1990) H I survey shows a cloud abutting the thickest part of the larger crescent. 
We were unable to find additional confirmation in radio synchrotron surveys (Haslam et al. 
1982; Alvarez et al. 97) or the R0SAT3/A keV soft X-ray survey (Snowden et al. 1997). The 
lack of confirmation is not a confirmation of the lack of a remnant for the following reasons. 
Caswell & Lerche (1979) empirically determined that, on average, the radio brightness of a 
SNR decreases with z. Correspondingly, our simulations of SNRs in low density environments 
found that their shells were not very dense and so probably not very emissive in synchrotron 
emission. In addition, the 3/4 keV emission from the simulated SNRs wanes far earher than 
does the 1/4 keV emission, thus a remnant can be observed in 1/4 keV X-rays after it has 
become too dim to be observed in 3/4 keV X-rays. The smaller feature, which is around 
I — 215°, b — —68°, also has an edge brightened crescent shape, though it is less crisp and 
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less bright than the larger feature. Suggesting that this is a SNR would be more speculative. 

The next element of study is the mottled background. Galactic X-ray emitting gas 
covers roughly half of sky above \b\ — 45° sky in Snowden et al. (1998)'s 1/4 keV survey map 
(assuming that the extragalactic sources contribute 400 xlO~^ counts s~^ arcmin"^). This 
compares well with the fraction of the sky covered by post shell formation phase SNRs (30 
to 90%). 

Finally, we turn to the spectral temperature of the emission. Several studies of ROSATs 
low resolution spectra, including Kuntz & Snowden (2000) 's and Snowden et al. (2000) 's 
found that the high latitude extraplanar emission has two temperature components, a ~ 
10^ K, non-homogeneously distributed component and a weaker, ~ 3x 10^ K, homogeneously 
distributed component. Because they are distributed differently, the 10^ K and the 3 x 10^ K 
components are thought to originate in different regions. If extraplanar SNRs are producing 
much of the 10^ K component, then we would expect their soft X-ray emission to look like 
a 10^ K spectrum. This was found to be the case when the simulated spectra were treated 
as real spectra (i.e. convolved with the ROSAT response function and fit to coUisional 
equilibrium spectral models), in Paper II. Note that the SNR plasma was often out of 
coUisional ionizational equilibrium and as a result, the gas's temperature was generally higher 
or lower than the temperature calculated from its spectrum. 



4.2. O"*"^ Column Densities and Intensities 

Soft X-rays are emitted by plasma whose carbon, nitrogen, and oxygen atoms have been 
ionized to the hydrogen-like and helium-like levels. The neighboring ions, lithium-like ions 
of carbon, nitrogen, and oxygen can be observed by absorption or emission in their strong 
resonance line doublets in the ultraviolet region of the electromagnetic spectrum (0+^: 1032, 
1038 A; N+4: 1239, 1243 A; 0+^: 1548, 1551 A). Owing to its large ionization potential, 
the observed O^^ is thought to be produced via coUisional ionization. N^"^, and C"*"^ can 
be produced via photoionization or coUisional ionization. In fact, Slavin et al. (2000) have 
already shown that substantial numbers of extraplanar ions have been ionized to C"*"^ 
by SNR photons. Our code simulates coUisional ionization and recombination but does not 
simulation photoionization. Therefore we will confine our comparisons to the soft X-ray and 
0+5 data. 

Hundreds of Milky Way sightlines have been surveyed for O"'"^ absorption (Jenkins & 
Meloy (1974); Jenkins (1978a); Hurwitz & Bowyer (1996); Widmann et al. (1998); Zsargo 
et al. (2003); Wakker et al. (2003); Oegerle et al. (2005)). The data suggest that this 
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ion is distributed inhomogeneously and is more plentiful nearer to the plane than at greater 
distances. The column density vs height data are plotted in the standard fashion in Figure 4. 
For comparison, our most productive model {Eo = 1.0 x 10^^ ergs, -Be// = 5/iG, drag 
coefficient = 0) and least productive model (Eo = 0.5 x 10^^ ergs, i?e// = 5fiG, nonbuoyant) 
have been plotted. Because our models begin at z = 130 pc, we have also plotted the 
expected local contribution (taken from Oegerle et al. (2005)) and added it to our models 
in order to obtain the total quantity of 0+^ ions expected to reside between the Sun and 
various heights above the plane. The extraplanar SNR column densities reasonably track the 
observations within the relevant height range (~ 130 < z <~ 2000 pc). Thus, extraplanar 
SNRs can explain both the observed column densities and the increase in column density 
with height, within this z range. 

The extraplanar SNR model was not meant to explain the 0+^ observed within the first 
hundred or so parsecs of the Sun. That material has been attributed to the Local Bubble 
and clouds within it. Nor does the model explain the O"*"^ residing beyond z ~ 2000 pc 
(about 6 times the SN scaleheight) . That column density requires another source. 
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SNR simulations (blue and purple curves) combined with the estimated Local Bubble con- 
tribution (Oegerle et al. (2005), dashed red line). The purple and blue curves correspond to 
the least and most prolific extraplanar SNR models, respectively. The squares mark Jenkins 
(1978a)'s Copernicus observations. The Oi^Fi??/!? observations by Hurwitz & Bowyer (1996) 
are marked with downward pointing triangles while the ORFEUS //observations by Wid- 
mann et al. (1998) are marked with upward pointing triangles. The FUSE Local Bubble 
observations by Oegerle et al. (2005) are marked with diamonds. The FUSE Milky Way 
halo observations by Zsargo et al. (2003) are marked with circles. The two FUSE Milky Way 
star observations by Wakker et al. (2003) are marked with stars. For the above plotting, 
detections are marked with open symbols and upper limits are marked by solid symbols. The 
Howk et al. (2002) SMC observations are marked with asterisks. The Howk et al. (2002) 
LMC observations of O^^ within the velocity range v < ±50 km s~^ are plotted with Xs. 
These column densities plus higher positive velocity O"^^ observed along the same sightlines 
are plotted with pluses. All of the Howk et al. (2002) observations were detections. Vertical 
lines connect the higher and lower LMC datapoints. The agreement between the model 
and the data suggests that the extraplanar supernova remnants reasonably explain the 0+^ 
observed between z ~ 130 pc and z ~ 2000 pc. 
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The upper limits and weak detections in Figure 4 demonstrate that the 0"'"^-rich gas 
docs not completely cover the high latitude, extraplanar sky. Of the sightlines terminating 
above 130 pc, 1/8 have O^^ upper limits and an additional 1/8 have such small column 
densities (less 0+^ than expected from local material) that we can assume they have not 
encountered extraplanar O"*"^. In comparison, our simulation results also imply incomplete 
(30 to 90%) sky coverage. 

Wakkcr ct al. (2003) 's Figure 9 displays the column density vs velocity distribution of 
the 0+^ features observed in their survey. The distribution consists of a cluster of features 
with —50 km sec~^ < f < 50 km sec^^ plus separate clusters having —400 km sec~^ v < 
— 100 km scc^^ and 75 km sec^^ ^ v < 400 km sec^^. The high velocity clusters are not 
symmetric; the positive high velocity cluster extends to lower speeds and column densities 
than does the negative high velocity cluster. In comparing with our simulation results, we 
must keep in mind that our predicted velocities are in the vertical direction while Wakker et 
al. (2003)'s velocities arc along high latitude lines of sight. Thus, we would want to multiply 
the observed velocities by cosecant (latititudc), where the typical latitude is around 50° to 
60°. In our calculations, the greatest velocity to which the extraplanar SNRs are accelerated 
is 140 km sec~^ and this is for only one set of parameter choices. Such remnants would 
contribute to the moderate speed, low column density portion of the observed 75 km sec~^ 
< V < 400 km sec~^ cluster, though they would not produce the entire cluster or its negative 
velocity counterpart. If extraplanar SNRs experience moderate to large drag, then their 0+^ 
features would contribute to the —50 km sec~^ < v < 50 km sec~^ cluster. While both 
scenarios lie within the broad range of possibility, the first scenario is highly constrained 
and thus unlikely. It is far more likely that the observed nearly symmetrical, low velocity 
distribution includes extraplanar SNRs who have experienced moderate to large drag. 

Now we turn to 0+^ resonance line emission. Observations of the resonance line intensity 
have the potential to nicely complement the column density data. The intensity produced by 
a column of 0+^ ions is sensitive to the local volume density and temperature. Emissivity 
goes as density squared. Hot gas emits readily, while warm/hot gas emits sparingly and 
cool gas does not emit. Using Voyager, Murthy et al. (2001) set upper limits on the O^^ 
doublet intensity for hundreds of sightlines. Their tightest was a 90% confidence upper limit 
of 2600 photons s"^ cm-^ sr'^ toward £ = 117°, b = 51°. Using FUSE, Shelton et al. (2001), 
Dixon, et al. (2001), Shelton (2002), and Otte et al. (2003) have detected 0+^ emission 
from a handful of high latitude directions outside of young SNRs and superbubbles. Their 
doublet intensities range from ~ 2400 to ~ 4700 photons s~^ cm~^ sr~^. Ignoring the effects 
of absorption and subtracting the Local Bubble contribution {2a upper limit = 800 photons 
s~^ cm~^ sr~^, Shelton (2003)) leaves 1600 to 4700 photons s~^ cm~^ sr~^ to be attributed 
to the extraplanar sky. In comparison, our model predicts 840 to 2100 photons s~^ cm~^ 
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sr~^, with the variation mostly due to assumed explosion energy. Given that the extraplanar 
SNRs could explain the O"*"^ column density between z = 130pc and 2000 pc, but not the 
higher O"*"^, the same could be true of the emission. To determine the answer, we will need 
to know the typical height of the observed 0+^ emission. 

5. Summary 

• This paper refines Papers I and II's preliminary calculations of isolated extraplanar 
SNRs by modeling the decrease in interstellar density with height above the Galactic 
plane, considering buoyancy, exploring a wider range of explosion energies and nonther- 
mal pressures and predicting a wider range of observables (O'*'^, N"*"^, C^^ emission and 
absorption column densities, 1/4 keV X-ray surface brightness, area coverage, volume 
occupation, vertical velocity, and variation with height above the Galactic plane). 

• The predicted average O"*"^ column density matches the observed distribution of column 
densities as a function of height above the Galactic plane. All 9 of our models (3 sets 
of explosion energies and ambient effective magnetic field strengths and 3 sets of drag 
coefficients) fell within the scatter in the observed A'"(0+^) sin |6| versus \z\ for the 
height range of interest {z up to 6 times the SN scaleheight). See Subsection 4.2. 

• Isolated extraplanar SNRs can explain the Galactic halo's 1/4 keV X-ray brightness, 
within the precision of the SN rate estimates. Assuming that these SNRs explode with 
energies of 1 x 10^^ ergs and experience an ambient effective magnetic field of 5 /xG, 
then the population of isolated SNRs born above the Galactic H I layer (scaleheight 
~ 130 pc) produces a time and space averaged 1/4 keV soft X-ray surface brightness 
of ~ 310 X 10-^ ROSAT counts s"^ arcmin-^. This is 80% of the observationally 
determined surface brightness of the Galaxy's southern high latitude sky beyond the H I 
layer (400 x 10~^ ROSAT coxmis s~^ arcmin"^). If we were to use a larger SN explosion 
rate, such as that estimated by McKee & Williams (1997), then the prediction would 
rise to 400 x 10^*^ i^O^ylT counts s^^ arcmin"^, accounting for the entire observationally 
determined surface brightness. However, using the rate adopted in this paper and the 
other choices of SN parameters yields average surface brightnesses as low as 120 x 10~^ 
i?(95'^T' counts s~^ arcmin"^. Buoyancy makes little difference here. See Subsection 
4.1. 

• The predicted average 0+^ intensity (840 to 2100 photons s~^ cm~^ sr~^) contributes 
significantly to the observationally determined intensity originating beyond the Local 
Bubble (~ 2400 to ~ 4700 photons s"^ cm'^ sr-^). See Subsection 4.2. 
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• Before SNRs develop a cool shell, they are bright in 1/4 keV emission, and 0+^, N+^, 
and C"*"^ resonance line emission, and are rich in O"*"^, N+^, and C"*"^ ions. After the 
shell forms, the remnants dim, but still harbor large numbers of lithium-like ions. 
Thus, the spatial distributions of 1/4 kcV emission and O"*"^ column densities should 
be de-correlated, as found by observers (ex: Savage et al. (2003) 's. Figure 9). 

• Old SNRs, which glow dimly in 1/4 keV X-ray and UV resonance line photons and 
retain observable numbers of lithium-like ions, should cover 30 to 90% of the high 
latitude sky. Young remnants, which glow brightly in these regimes should cover less 
than a percent of the high latitude sky. Correspondingly, maps of the observed soft 
X-ray background beyond the Galactic H I layer, outside of supcrbubblcs, and after 
subtracting the extragalactic component appear mottled, as if roughly half the high 
latitude sky is covered with dim emitting regions while the other half is bare. The 0+^ 
column density surveys also find null detections and low detections. See Subsections 
3.5, 4.1, and 4.2. 

• The fraction of volume occupied by hot supernova remnant bubbles peaks around the 
supernova scaleheight and is less than 10% at all heights. The volume filled by young 
remnants is only a fraction of a percent. See subsection 3.2. 

• We calculated the vertical displacement and velocity distribution, assuming buoyant 
acceleration and zero to moderate drag. In the extreme case of unfettered buoyancy, 
the velocity can reach -1-140 km sec~^, but if SNRs experience modest drag, then their 
velocities will range from to 30 km/sec. In order for the gas to rise and fall again, 
the drag coefficient must exceed 1.0. Buoyancy is not absolutely necessary to explain 
the observed vertical distribution of 0+^ ions. See Figure 1 and Section 3.3. 

• Of the explored range of parameters, we prefer Case 3 {Eo = 1.0 x 10^^ ergs, -Be// = 
5.0/iG) because its predictions best match the high ion and soft X-ray data. The 
predictions from the other cases are factors of ~ 2 to ~ 3 lower. Of the explored range 
of drag coefficients, we prefer moderate to large drag, because it limits the buoyant 
velocity to 50 km sec~^ or less. High latitude observations find many O"*"^ features in 
this velocity range. 

• This picture can be tested by: 1.) examining the high latitude crescent shaped bright 

X-ray regions (/ ~ 247°, b 64° and I ~ 215°, b 68°) to determine if they are the 

hmbs of young remnants (see Subsection 4.1), and 2.) and testing the spectral features 
of the mottled dim 1/4 keV background for signatures of old, cooling, recombining 
SNR gas. (See Paper II.) 
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